Critical bending point in the Lyapunov localization spectra of many-particle systems 
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The localization spectra of Lyapunov vectors in many-particle systems at low density exhibit 
a characteristic bending behavior. It is shown that this behavior is due to a restriction on the 
maximum number of the most localized Lyapunov vectors determined by the system configuration 
and mutual orthogonality. For a quasi-one-dimensional system this leads to a predicted bending 
point at n c w 0.432iV for an N particle system. Numerical evidence is presented that confirms this 
predicted bending point as a function of the number of particles N. 
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The Lyapunov spectrum is an indicator of dynamical 
instability in the phase space of many-particle systems. It 
is introduced as the sorted set {X^} n , A^ > A^ 2 ^ > ■ • • 
of Lyapunov exponents A*-™-*, which give the exponen- 
tial rates of expansion or contraction of the distance 
between nearby trajectories (Lyapunov vector) and is 
defined for each independent component of the phase 
space. In Hamiltonian systems and some thermostated 
systems, a symmetric structure of the Lyapunov spec- 
tra-the so called the conjugate pairing rule, is observed 
[H El- B- One of the most significant points of the 
Lyapunov spectrum is that each Lyapunov exponent in- 
dicates a time scale given by the inverse of the Lya- 
punov exponent so we can consider the Lyapunov spec- 
trum as a spectrum of time-scales. The smallest posi- 
tive Lyapunov exponent region of the spectrum is dom- 
inated by macroscopic time and length scale behavior, 
and here some delocalized mode-like structures (the Lya- 
punov modes) have been observed in the Lyapunov vec- 
tors [!IE0HH[i![ilHm0. 0n the other hand ' 

the largest Lyapunov exponent region of Lyapunov spec- 
trum is dominated by short time scale behavior, and in 
this region the Lyapunov vectors are localized (Lyapunov 
localization) . The position of the localized region of Lya- 
punov vectors moves as a function of time. A variety of 
many-particle systems show Lyapunov localization, for 
example, the Kuramoto-Sivashinsky model 1 1511, a, ran- 
dom matrix model 0], map systems [l7t llct Il9l l20| . 
coupled nonlinear oscillators 21], and many-disk systems 

mm. 

Recently, a quantity to measure strength of Lyapunov 
localization was proposed [2^ . For each Lyapunov expo- 
nent we construct the normalized Lyapunov vector com- 
ponent amplitude 7^ = |<5r^™' 1 |/|<5r^™- ) | for each particle 
(™) rr(") 



Here ST (n) = (Sr[ n \ST 



, Srffi ) is the Lya- 
punov vector for the n-th Lyapunov exponent \( n \ and 
ST^ is the contribution of the j-th particle to the n-th 
Lyapunov vector. The localization of the n-th Lyapunov 
vector is then 
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The bracket (• • • ) in Eq. indicates the time-average. 

The quantity — Y^Liilj m 7j™'} m the definition (JIJ of 

can be regarded as an entropy-like quantity, as 7^™ 
is a distribution function over the particle index j. The 
quantity W (n) satisfies the inequality 1 < W (n) < N, 
and can be interpreted as the effective number of particles 
contributing to the non-zero components of the Lyapunov 
vector. In a Hamiltonian system it satisfies the conjugate 
relation y\A- D_J+1 ) = VVW') for any j with the phase space 
dimension T>, because of the symplectic structure. 

The set of quantities {W^} n which we call the Lya- 
punov localization spectrum, has been calculated previ- 
ously in many-particle systems with hard-core interac- 
tions [2^] and with soft-core interactions |l3|. In these 
systems, the value of usually increases with the Lya- 
punov index n, and this implies that Lyapunov vectors 
for the largest exponent region are the most localized. 
The quantity W (rl) has a minimum value of 2, as a min- 
imum of two particles are involved in each collision, and 
it has been shown numerically that the value of yy(") for 
the largest Lyapunov exponent approaches 2 as the den- 
sity approaches zero (23J. It was also shown that W'™- 1 
can detect not only the localized behavior of Lyapunov 
vectors, but also the de-localized behavior observed in 
the Lyapunov modes p^ ]. 

One of the important characteristics of the Lyapunov 
localization spectrum is its bending behavior at low den- 
sity [2i|]. In Fig. [I]we show an example of such bend- 
ing behavior |24| in a quasi-one-dimensional system of 50 
hard disks with periodic boundary conditions. Here the 
radius of the particles R = 1, the mass of each particle 
m = 1, the total energy of system E = N, and the system 
size L y = 2i?(l + 10~ 6 ) and L x = NL y (l+d) with param- 
eter d, related to the density by p = nR 2 /[(l + d)L y ]. In 
this system the particles remain ordered because the ver- 
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tical size L y prohibits the exchange of particles [8|, .12, 2 
In Fig. ^ five different Lyapunov localization spectra 
are represented at different densities. In the low density 
limit, this figure shows a bending point in the Lyapunov 
localization spectra at n c . At this bending point n w n c 
the Lyapunov localization spectra changes from a 

linear dependence upon the Lyapunov index (n) to an ex- 
ponential dependence. This bending becomes sharper at 
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FIG. 1: (Color online) The normalized Lyapunov localiza- 
tion spectra {(W' n - ) — 2)/N} n as functions of the normalized 
Lyapunov index n/(2N) for quasi-one-dimensional 50 hard- 
disk systems at different densities given by d — 1 (diamonds) , 
d = 10 (circles), d — 100 (triangles), d — 500 (inverted trian- 
gles), and d = 10000 (squares). The parameter d is inversely 
proportional to the density as p — irR 2 /[(l + d)Ly]. The solid 
line is the asymptotic form of Lyapunov localization spectrum 
as the density approaches zero, and n ~ n c is its bending 
point. 

lower density, and the numerical results suggest that in 
the low density limit the Lyapunov localization spectrum 
converges to the solid line given by 

w („) _ f 2 in n < n c , . 

\ -y-aNexp{-f3n/(2N)} in n > n c y > 

with constants a, f3, and 7 = 2 + aN cxp{— f3n c /(2N)} 
[so that W (n) | n= „ c = 2 in Eq. ©]. In other words, 
we can estimate the critical value n k, n c defined as the 
bending point of the Lyapunov localization spectrum by 
the value of the fitting parameter n c in the fit of the 
Lyapunov localization spectrum to the function J3J. 

The bending behavior of Lyapunov localization spec- 
tra, shown in Fig. ^is associated with a similar bending 
point in the Lyapunov exponent spectra. Moreover, the 
existence of a linear dependence of yy( n ) on n in the re- 
gion n < n c , at low density, is connected with some other 
known kinetic properties, for example, that the mean free 
time is inversely proportional to density, and the Krylov 
relation, that the largest Lyapunov exponent A^ 1 ) de- 
pends on the density p like \W ~ —p\np [2j|. These re- 
sults suggest that the existence of the linear dependence 
of is connected to the density range where kinetic 

theory provides an accurate description. These points 
were investigated in detail in quasi-one-dimcnsional sys- 
tems, although the bending behavior of Lyapunov local- 
ization spectra is also observed in fully two-dimensional 
square systems |23f . However, no mechanism was pro- 
posed for this bending behavior. 

In this paper we construct a mechanism that leads to 



the bending behavior observed in the Lyapunov localiza- 
tion spectra, and predicts the critical value n c for quasi- 
one-dimensional systems. 

We begin by recalling some properties of the Lyapunov 
vectors of many-particle systems. The first property ex- 
presses the fact that different Lyapunov vectors sample 
different directions in phase space. 

[I] Lyapunov vectors with different Lyapunov indexes are 
orthogonal. 

The second property of Lyapunov vectors is based on the 
fact that particle interactions occur between two different 
particles: 

[HA] In the low density limit, all the Lyapunov vectors 
SF^ n \ for n < n c have non-zero components for 
only two particles. 

This leads to the known result that — > 2 as p — > 

0. The third property is justified only for quasi-one- 
dimensional systems, and restricts the property [HA] fur- 
ther to 

[IIB] Two particles, whose Lyapunov vector components 
of ST^ n ', n < n c are non-zero values in the low 
density limit, are nearest-neighbors. 

In Ref. [2i| it is shown that in the largest Lyapunov 
exponent region, at low density, two non-zero particle 
components of the Lyapunov vector appear by particle 
collisions. Moreover, in quasi-one-dimensional systems, 
particle collisions occur only between nearest neighbor 
particles. Based on these facts, the property [IB] is jus- 
tified. 

Using the properties [HA] and [IB], in the low den- 
sity limit, the Lyapunov vectors 6T^ n ' in the largest Lya- 
punov exponent region can be represented as 

$r™ „ ( o, 0, • , 0, jr<^ar£ } +1 ,o, • • • , o) (3) 

where is the null vector. Here, the particle number- 
ing for the quasi-one-dimensional system is from left to 
right, and fi n and fi n + 1 are nearest neighbor particles 
whose Lyapunov vector components are non-zero. (We 
put ST^ +1 = dT^ for periodic boundary conditions in 
the longitudinal direction.) In general, the particle num- 
ber p n depends on the Lyapunov index n and on time. 
Now we consider the restriction imposed by condition [I] 
on the form (|3J| of Lyapunov vector. To satisfy condition 
[I], the particle numbers p n and fi n + 1 in the Lyapunov 
vector (0 must be different for different Lyapunov in- 
dexes in n < n c . Here, we assume that the vector com- 
ponents for particle j, non-zero components ST^ and 

<5r^ n ' are themselves not orthogonal for different Lya- 
punov indices n and n'. This restricts the number of 
independent Lyapunov vectors of the form J3J), and puts 
an upper limit on the critical value n c . 

This situation is explained using a simple model whose 
schematic illustration is given in Fig. [3 In this randomly 
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FIG. 2: A schematic illustration of the randomly distributed 
brick model, used to explain the maximum number of in- 
dependent Lyapunov vectors with only two non-zero particle 
components. The box labelled j corresponds to the Lyapunov 
vector component of the j-th particle (j — 1, 2, • • • , N). Gray- 
filled rectangular bricks occupying boxes fj, n and fi n + 1 rep- 
resent possible non-zero components of Lyapunov vectors for 
which n < n c . The value n c is the average number of bricks 
that can be randomly arranged on the set of boxes. 
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FIG. 3: The normalized critical value n c / (2N) given by the 
randomly distributed brick model as a function of the number 
of particles TV. The solid line is given by fitting the numerical 
data to a constant function n c / (2iV) = £ with the fitting 
parameter £ w 0.216. 



distributed brick model, N boxes corresponding to each 
of the particles are arranged on a line. Each Lyapunov 
vector represented by Eq. © must have non-zero com- 
ponents for only two particles ST^ and 5T^ +1 . These 
are shown as a gray-filled rectangular brick filling boxes 
fx n and (fi n + 1) in Fig. |21 To satisfy condition [I], these 
bricks must not overlap. The critical value n c is then 
the average number of bricks which can be put without 
overlaps on the /V different boxes. 

There is one more important point to get an explicit 
value of n r from the above mechanism: 



[HI] The particle number ji n in Eq. Q is chosen ran- 
domly with respect to the Lyapunov index n(< n c ), 
so that there is no overlap among the non-zero 
Lyapunov vector components for Lyapunov vectors 
with different Lyapunov indices. 

This means that the bricks shown in Fig. [21 must be 
randomly distributed without overlaps. Therefore ran- 
domly constructed configurations with any number of 





0.5 




0.45 




0.4 




0.35 








0.3 


c\T 

i 


0.25 




0.2 




0.15 




0.1 




0.05 








no/(2N) 



f A 
N = 40 ♦ 
N = 50 - 
N = 60 ▼ 
N = 70 B 

V J 



0.1 



0.2 



0.3 



0.4 



0.5 



0.6 



n / (2N) 



FIG. 4: (Color online) The normalized Lyapunov localization 
spectra {(W'™' — 2)/N}„ as functions of the normalized Lya- 
punov index n/(27V) in quasi-one-dimensional systems of dif- 
ferent size; N = 40 (diamonds), N = 50 (triangles), TV = 60 
(inverted triangles), and N = 70 (squares) for a density given 
by d = 10000. The lines are the asymptotic forms of Lya- 
punov localization spectra in the low density limit for these 
different values of N. 



single empty boxes at non-neighboring positions are pos- 
sible. Taking an ensemble average of the values of n c for 
each possible configuration we obtain the critical value 
n c . Figure 13 shows the normalized critical value n c /(2N) 
from such a numerical simulation for different numbers of 
particles N. The data in this figure shows that n c /(2N) 
is independent of N, giving an estimate of the normalized 
critical value as 



2N 



0.216. 



(4) 



The accuracy of this value can be used as a check of the 
proposed mechanism for the bending behavior of Lya- 
punov localization spectra at low density. 

To draw the solid line shown in Fig. \\\we have already 
used the value obtained in Eq. The parameters a 

and f3 in Eq. J2J are regarded as fitting parameters, and 
we used values of (a, (3) ~ (2.51, 7.26) for the line in Fig. 
n The result shown in Fig. ^ using the critical value Q 
gives a very satisfactory fit. 

In order to check that the result (0J is satisfied for 
any number of particles N, we present Fig. It is 
the [n/(2iV)]-dependence of [W (n) - 2}/N in quasi-one- 
dimensional systems of different sizes at the same density. 
Notice that values of [W'™' — 2]/N themselves decrease 
slightly as the number of particles increases, meaning 
that the fitting parameters a and (3 in Eq. @ depend 
only slightly on N. The values of the fitting parameters 
a and (3 for Fig. H were (a, (3) w (2.83, 7.83) for N = 40 
(dotted broken line), (a,/3) « (2.51,7.26) for N = 50 
(dotted line), (a,/?) « (2.28,6.80) for N = 60 (broken 
line), and {a, (3) w (2.10,6.39) for N = 70 (solid line). 
However, all the data for different numbers of particles 
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in Fig. 0] are nicely fitted using the same critical value 
n c /(2N), given by Eq. (QJ. 

To conclude, we have shown that a model of randomly 
distributed bricks on a line (Fig. [2J can predict the max- 
imum number of Lyapunov vectors which have a Lya- 
punov localization equal to two, in the low density limit. 
These are the Lyapunov vectors for which n < n c in the 
linear region that lead to the bending behavior of Lya- 
punov localization spectra. We showed that this behavior 
comes from a restriction on the maximum number of the 
most localized Lyapunov vectors with non-zero compo- 
nents for only two particles. The randomly distributed 
brick model was applied to quasi-one-dimensional sys- 
tems, giving the specific value n c /(2N) ks 0.216 for the 
critical value n c of the Lyapunov localization spectra for 
any number of particles N. Numerical simulations give a 
value of n c /{2N) w 0.22 ± 0.06. Our explanation for the 
bending behavior is independent of the system width L y , 
as long as the particle order is invariant and the sys- 



tem remains quasi-one-dimensional. We have checked 
this using L y = 2R(2 — 10~ 6 ) and the critical value is 
unchanged. 

The critical value of Lyapunov localization spectra de- 
pends upon the shape of the system. Numerical re- 
sults show that the critical value of the Lyapunov lo- 
calization spectra for a square system is smaller than 
that for the quasi-one-dimensional system 23] . The ran- 
domly distributed brick model is specific to the quasi- 
one-dimensional system and would need to be generalized 
for a square system, and property [EB] can no longer be 
assumed. The calculation of the critical value n c for a 
general two (or three) dimensional systems remains as a 
future problem. 
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In the normalized form /N}„ of Lyapunov local- 

ization spectrum, its minimum value 2/N depends on 
N. To effectively compare Lyapunov localization spectra 
with different numbers of particles, we use the normal- 
ized Lyapunov localization spectra minus 2/N in Figs. 
and |2| This is crucial in Fig. [5] where the particle number 
dependence of Lyapunov localization spectra is discussed. 



